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ABSTRACT 

Two approximate interval estimation procedures for mechanical component reli- 
ability, P(X > Y), are-developed and their accuracy evaluated by computer simulations. 
The strength, X, of the component and the stress, Y, applied to it are independent 
normally distributed variables with unknown means and variances. In-the first interval 
procedure the variances are equa!. In the second procedure the variances may be une- 
qual. 

The derived intervals are quite accurate for the cases simulated which include large 
and small sampis sizes. These procedures are simple to apply and require the use of 
percentile points of the Student’s t distribution. In the second procedure, the degrees 
of freedom of the associated t statistic is a function of the test data , and therefore it is 


random. 
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I. INTRODUCTION 


Let X and Y be independent random variables with normal cumulative distribution 
functions F,(x) and F,{y)-respectively. Suppose X is the strength of a mechanical com- 
ponent, and Y is the stress applied to the component. The strength depends on the 
material properties, manufacturing, procedures.and other factors. The stress is a function 
of the environment to which the component is subjected. Component failure-is defined 
by the event Y < }'. Component reliability,-R, is defined by R = P(X > ¥). Ris called 
mechanical reliability. Two lower-confidence limit procedures for R are developed in this 
thesis. In both procedures the means and variances of X and Y are unknown. In one 
procedure the variances are-assumed to be equal. 

A nonpararvtric interval estimation procedure for R was first proposed by 
Birnbaum [Ref. 1: pp. 13-17] using the Mann-Whitney U statistic. Birnbaum and 
McCarty developed a procedure for computing the minimum sample size needed for 
such a confidence imierval to have a given width and confidence Jevel [Ref 2: 
pp. 558-562]. Owen, Craswell and Hanson (Ref. 3: pp. 906-924] provided more detailed 
tables for use in computing sample sizes and confidence intervals for the 
Birnbaum-McCarts procedure. Tables designed especially for the normal distribution 
were also included in their paper. Govindarajulu (Ref. 4: pp. 229-238] observed that the 
bounds employed by Birnbaum and McCarty for obtaining the confidence intervals can 
be substantially improved asymptotically and reduced the Birbaum-McCarty bounds by 
approximately 1,2. Church and Harris (Ref. 5: pp. 49-54] pointed out that the sample 


sizes required by these nonparametric procedure are likely to be too large for many 


practical situations. 











Owen [Ref. 6: pp.445-478] gave an exact confidence limit procedure for 
R= P(X =x) where X is normally distributed with unknown mean and unknown vari- 
ance, and x is a-constant. His-procedure uses the noncentral t distribution and extensive 
table-lookups are needed. Owen-and JJua [Ref. 7: pp. 285-311] developed special tables 
that reduced these-calculations. Their tables were limited to two confidence Jevel values; 
namely 90° and.95%. Lee [Ref. 8: pp. 15-22] reports on a closed form equation for this 
confidence limit that is approximate but quite accurate. His equation applies for any 
confidence level and-uses the central student’s t distribution. 

Church and Harris [Ref. 5: pp. 49-54] developed approximate confidence limits for 
R= P(X> Y) , under the assumption that the stress, Y, has.a standard normal distrib- 
ution and F,(x) is normally distributed with unknown mean and variance. 


Throughout this thesis we write Y ~ (nu, o?) to denote that X has.a normal distrib- 


ution with mean uw and variance o?. 








Hi. APPROXIMATE INTERVAL ESTIMATION PROCEDURE FOR 
RELIABILITY R = P(A > ¥) — EQUAL VARIANCE CASE 


A. LOWER CONFIDENCE LIMIT PROCEDURE 

Suppose the strength, \, of a mechanical device and the stress, Y, applied to it are 
‘independent variables, with normal probability distributions. We assume both means 
are unknown and buth variances are unknown but have common values o?. The me- 


chanical reliability, R, of the device is defined as follows: 














R=PiX>¥J 
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where @ is the standard normai cumulative distribution function. Let 6=——3—. 
Then 
R= tS}, (2.2) 
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A consistent estimator [Ref. 9: PP.289-294] of 6 is | 
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is the pooled sample variance of X and Y; n is the size of the sample on X; mis the:size 
of the sample on- Y;-and ¥ and Y are the respective sample means. Since X and Y are 


independently normally distributed, . 


¥-¥~ Nuy-ayo( t+ ) : 


The general method for deriving confidence intervals [Ref. 9: PP.347-355] can be 


‘used to obtain a lower 100(1 — «)% confidence limit, ee , for 6. Suppose as 





ar. 
amis 
constructed. from the data. Then under the general method, 6,.,_, is the value of 6-such 
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where N.C.T. denotes a noncentral t random variable with noncentrality parameter 


6 Myo By Xx~-y 
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Jf free Afin 





The degrees of freedom of this noncentral t is n+m-2. One can use the noncentral t 


table to obtain the solution for 6 in Equation (2.4) which is bn, j-2 + Owen and Hua [Ref 








7: pp. 285-311] have taken this approach for the univariate case. Their method-required 
the extensive use of tables. We shall take a different approach ‘here, in order to find a 
closed expression for an approximation to b mn 

In his thesis, Lee (Ref. 8: pp. 15-22] developed an approximate 100(1 — «)%o lower 
confidence limit for R= P(Y > x) where X ~ N(u, 0”), w and o? are unknown, and x a 


constant. Lee’s expression is 





2 
A » td Ky 12 : 
67 peg = Kym Ny, ep 2:5 
L, l=x x E is | Ina, n~] (2-5) 
where K, = A > ~ and tes, iS the 100(1 — «)" percentile point of the student's t dis- 


tribution with n-1 degrees of freedom. 

A method analogous to Lee’s procedure can be used to develop an equation for 
Bris in the bivariate case described at the outset of this section. If we substitute 
(1/n + 1/mi) for (1/2) and (+ a — 1) for (x)-in Equation (2.5), we would obtain the lower 
confidence limit 


2 
5 =k BERS a pt eects (0 1/2 
OL, Jaa = K— ( nt )+ n+ m—1—Jk) | U-x, n-m=2" (2.6) 


where K= 25+ The corresponding 100(1—a)%o lower confidence limit for 


R=PY> ¥)is 


A 
a 6p l~x 
=o, ——— |], 2.7 
Rr, l=x 4 ( a/2 } ( ) 


Our computer simulation results show that this lower confidence limit is quite accurate. 
The results are tabulated in Tables 1 and 2. A description of the simulation procedure 
along with a analysis on computer results will be given later in this chapter. The im- 


portant point of these tables is the comparison between R and Rigas, za-2 Lf these two 














values are equal, the approximate lower confidence limit procedures given in Equations 


(2.6).and (2.7) are nearly exact. 














Table 1. ANALYSIS OF 90% CONFIDENCE LIMIT APPROXIMATION OF 
, R = P(X> Y) FOR EQUAL VARIANCES CASE WITH EQN. 2.6 


Mean error Variance of 
— from R error 
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The above method for deriving ,,,., canno. be used to find a lower confidence in- 
terval for P(X > ¥)avhen we drop the assumptivn of equal variances. Consequently, 


we shall-develop 6 using a different approach which has greater potential for con- 


hey bon 


-‘structing-confidence.intervals when variances are ‘not equal. Let 








K=g(¥-¥, S)=tehS As ae (2:8) 
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The Taylor expansion of K at (uy— Hy, 0?) , using only first order derivatives, is given 


by 
T oT ef _ pase aya TY A F lye ny 08s 3 
aX — ¥, S*) = g(uy— ny 6°) + LX — ¥ - (uy xy] a =F) | een ane 
sigh eRe ; 
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The expected value and variance of g(\ — ¥, S?) are as follows: 


uy a by 


E(A) =4,= 3 
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An estimator for cy, is then 


(4 ty er 
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where K = $ 


Com py 
——— js close tothe 


For large sample-size, n and-m, the probability distribution of 
Ox 


‘standard normal distribution. An.approximate 100(] — «)%o lower confidence limit for 
Uy is K—G,Z,_,, where Z,_, is the 100(1 —-c)®0 percentile of the standard normal dis- 
: : : : p sch ; K- Uy ‘ Pane? F 
tribution. We choose to approximate the distribution of --,—— with a . distribution: 

6. 
that has n+m-2 degrees of freedor. This approximation should accommodate sinall 
samples better than the normal approximation. The computer simulations will reveal 


the-accuracy of this choice. Consequently an-approximate lower confide:.ce limit for 


Hy is-given by 


A A 
HRY ia K- CRY ~y, n=m-2 


=x-[(445 aa? +m=— 2) Sex, n=m—2 (2.12) 
A 
= Op, jx 


The corresponding 100(1 — a) % lowe.’ confidence limit, R,,.,, of component reli- 


ability is 


1. Example 
We illustrate the application of this ; *-edure with an example. We cempute 
the 90% lower confidence limit of the consponent reliability, R,.95 given the following 


data: 


Se 


or 


a 


a 





X: 9.26 10.19 9.79 ‘14.27 10.06 9..22 9.75 8.46 9.79 9.52 


11.55 9.41 9.99 9,22 11.22 9.89 9.54 9.24 8.77 9.86 








is 10.03 10..33 10.95 8.33 10.85 


¥: 8.07 8.53 7.74 325 32.5% 9.36 6.82 8.65 7.00 7.65 


7.94 9.15 6.99 7.Jc 5.88 9.56 9.58 7.85 8.91 7.52 


9.97 8.76 7.01 9.9€ 7 45 9.05 9.54 10.72 7.86 7.87 


m= 30, ¥ = 8.308 
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Oa. near ar = 3.978 


ubm—2 





Be cat 1.4 whi 
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A 
P Sia P 
Re oo= o( 102 = 0.798 
7 f2 


We note that the point estimator of R, 








a( 252) =of _ ) = (1.122) = 0.866, 


and the computed 90° lower confidence limit of R, R,.99 = 0.798 <.0.868 as it should 
be. We cannot draw any conclusion about the accuracy of any confidence interval 
procedure from one example. We need computer simulations to do tits. The next sec- 


tion considers the accuracy of this procedure. | 











wn 














B. ACCURACY CF THE PROCEDURE 
1. ‘Measvres of accuracy and the concept of.compurer simu!2tion 
The accuracy of the interval proceaure in Equation (2.12) was evaluated in 


terms of the following four characiéristics: 


¢ The actual confidence level of the interval, i.e. the portion of times an estimated 
lirrit will: cover the truc reliability. 


¢ The mean error between th: -estimated limit and the true reliability, which is de- 
noted as ‘mean error from R in the simulation result tables. 


« The variance of the error between the estimated limit a..d true reliability, which is 
-denoted as ‘variance of error’ in the tables. 


¢ The 100(1 — a)" percentiie point of the distribution of R, ,.,. 


The actual confidence level of an approximate confidence interval can be ac- 
cessed during a computer-simulation which we will discuss later. 
To campute the 100(1 — a)" percentile point of the R,.,.,, one-only need to ex- 


amine the definition of Ris as a lower 100(1 — «) % confidence limit for R, i.e. 


PR g<R) 1a (2.14) 


This equation says that, R is the 1G0(1 -- a) percentile point of the probabii:ty distrib- 
ution of R,, i-w = hus, if we construct the distribution of Reins by computer simulation, 
we should find that the 100(1 — «)* percentile point of our constructed distribution is 
R, provided R, js iS a true 100(1 — «)% lower confidence limit for R. Let Rj... denote 
the 100(1 — a) percentile poin: of Rees Then the quantity |; ,_,—- R] is a measure 
of the accuracy of the procedure. 

We can construct the distribution of Rs. by generating a large number, say 
1000, of random observations on Rian for a given set of parameter values 
, M1, By By, and R. The 100(1 — a) empirical percentile point of the distribution of 


Ry, is the 1000(1 — a)" ordered statistic of Re feat Rea Ltles} ° 





2. Computer-simulation 
a. Simulation procedure 
; Eighteen sets of values of 1, m, wy, Hy, and o? were chosen to perform the 
simulations. They were selected in a manner so that R= o( £52" | , for R = 0.90, 
Cc 
0.95, 0.99. Thus when random samples of X and Y are peieitted: the reliability 


R= P(X> Y) will be at designed v."1e. Parameters were also chosen according to the 


following rules in order to cover practical conditions: 
© Ux> Ky 
* Reliability, R- :-0.90, 0.95, 0.99 
e Standard Deviations, ¢: 1.0, 20.0 
* Sample sizes,-(n, m): /8, 8); (8, 30); (20, 30) 


Actual sets of parameters ar - tabulated on Appendix B. 

Each set of parameters describes a case. For each-case, the set of parameters 
. were used to-generate samples of size n and_m for-normal variates X and Y respectively. 
The developed methods were used.to estimate the-confidence limit for confidence levels 
of « = 0.05, 0:1, and 0.20. For each. individual case, this procedure is replicated 1000 
times producing 1000 random-observations of R,. i» ©Lhe L000(1 — «) ordered statistic 
of R,, j-. iS compared with the true reliability. The actua: confidence levels is computed 
‘by counting the number of the 1000 Ry ns statistics that fall below the true reliability. 
We.also computed:the sample mean and variance of the estimation error using the 1000 
generated values of Riis 


The simulation procedure can-be summarized as follows: 


1) Generate n random normal variates X, Y ~ N(u,o7); m random normal variates 
Y, ¥~ Muze’). 


2) Compute X, Y, S. 


3) Compute K, é ice 


A 


P Dic 3 
4) Compute KR, )., = of a fora = 0.05, 0.10, 0.20. 


ala 


13 








5) Repeat steps 1 through 4, for 1000 times. Then order the R,,., to get 
Rarta-ay Reenen s+ Riooou-ataesy » 1 Raoooy.za-a» from the smallest to the largest. 


6) Print Ryeooa-s), i102" 
7) Print max{Ro, 0-2 Ro.za-y SR} the actual-confidence level of the approximation 


. 7 t 
will be To00" 


§) Compute .and print the sample mean and variance of the estimation error, to 
measure the precision and stability of the approximation. 


b. Simulation language 
The cimulations were conducted on the N.P.S. mainframe IBM 3033 com- 
puter. The programming language used for the simulations is VS FORTRAN 2. The 
random number generator LLRANII was used-to generate normal variates. The IMSL 
statistics function TIN was used to compute the percentile of the student t distribution 
and ANORDF was used to compute the probability of standard normal distribution. 
The FORTRAN source code for the simulation on the approximate procedure are at- 
tached in Appendix A. The FORTRAN code for the simultaneous comparison simu- 
Jation of the approximate procedure and the nonparametric procedure are attached in 
Appendix C. 
3. Analysis of simulation results 
The simulation results are tabulated in Tables 3 through 5. The results shows 
that Rieke: ii-ay 1S Very close to the truce reliability for every case of simulation, so that 


we have developed a ver} exact procedure. The procedure is more nearly exact for large 


sample sizes. 











Table 3. sy be OF 80% CONFIDENCE LIMIT APPROXIMATION OF 
= P(\> Y) FOR EQUAL VARIANCES CASE WITH EQN. 2.12 


ee p True confi- Mean error Variance of 
; 1001-2), Lulz) dence level | from R error 
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Table 4. ANALYSIS OF 90% CONFIDENCE LIMIT APPROXIMATION OF 
R= P(X> Y) FOR EQUAL VARIANCES CASE WITH EQN. 2.12 


: an True confi- Mean error | Variance of 

fee || | eis | eee | ome | are) 
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Table 5. ANALYSIS OF 95% CONFIDENCE LIMIT APPROXIMATION. OF 
Re PIX? Y) FOR EQUAL VARIANCES CASE WITH EQN. 2.12 _ 


Mean error | Variance of 

Ryosoa—s9. t01-29 dence level from R 
1942 | 0128 
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Ill. NONPARAMETRIC DISTRIBUTION FREE LOWER CONFIDENCE 
BOUND PROCEDURE 


The distribution.free confidence bound procedure suggested by Govindarajulu [Ref- 
4: pp. 229-238], may be applied to our problem. However we should note that the 
Govindarajulu procedure requires large sample ‘sizes and the true reliability values re- 
moved from 0 and I. For comparative purposes only, we evaluated his-procedure using 
the same computer simulation and analysis methods that were performed on our devcl- 
oped procedure. The same data was used to evaluate both procedures with emphases 


on ime ta-a) and ‘mean error’. The results-are displayed in Tables 6 and 7. 




















Table 6 COMPAF ‘DN OF DIFFERENT PROCEDURES ON 90% LOWER 
2 CONFIDE 7%é LIMIT ESTIMATION OF R = P(X>Y), EQUAL 


VARIANCES CASE. 
_e oe ee Bound 


ES) oximate Estimation 
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Table 7. COMPARISON OF DIFFERENT PROCEDURES. ON 95% LOWER 
CONFIDENCE LIMIT ESTIMATION OF R = P(X>Y), EQUAL 


VARIANCES CASE. 
Approximate Eotimation _Nonparamett Bound 


[eT e [ a ae rte ne 
1000t—2), LUl-2) Pass Li=2) 

is Fatal — 8989 0866 | 7955 fr 65a 
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ilo see ieee 
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' IV. APPROXIMATE INTERVAL ESTIMATION PROCEDURE FOR 
RELIABILITY R = P(X > Y) — UNEQUAL VARIANCES CASE 


A. LOWER CONFIDENCE LIMIT PROCEDURE 
Let X denote-component strength where X ~ A(uyo2) Let Y denote stress applied 


to the component where Y ~ N(zyo?) . Then 
7 ta y, 2 2 
X— Y~ Npy— By, oy + oy). 


The component reliability is defined as follows: 


R=PLX>¥] 
=p Ap ¥= (ux = 4y) os ER 
Jot oy al oh +o} (4.1) 
= & HS a ) 
NAS eae 


where © is the standard normal cumulative distribution function. Let 





_ ey Hy 





(4.2) | 
ae Re 
Vox roy 
then R = (6). 
A consistent estimator of 6 is 
yey 
po eT a (4.3) 


ig 2 
V Sit S¥ 


where S3 and S3 denotes sample variances and Y and Y are the respective sample 


means. 


2] 











xP 
For.an observed value of —==—— 


Sh sh 
6, using the general confidence interval method, is the value of 6 in Equation (4.2) such 


, the lower 100(1 — a) % confidence limit b, for 


that 


peep ate 


<——=— |. (4.4) 


The probability statement in (4.4)-cannot be reduced to an equivalent statement 
about a random variable whose distribution has been tabulated. This problem is similar 
to the Behrens-Fisher problem of finding a confidence interval for n,— py When both 
variances o} and o} are unknown and unequal. B. Lk. Welch (Ref. 10: pp. 28-35] has 
proposed .an approximate confidence interval for ny —-wy in this case using the statistic 
= . Welch approximated the distribution of this statistic with a Student's t 


distribution with v degreés of freedom. The degrees of-freedom chosen by Welch is given 


by 


° 


7 (S3fa + SPfmiy? 
S¥in? Sha” 


n— m=] 


The Welch statistic.is not useful in our problem, because we need to use a statistic that 


is a consistent estimator for —t=—24- 
¢ ¢ rarer a . 
Vox toy ay, eS 
In particular we choose to use the statistic g= se The mean and variance 
5% + 5} 
F en : ah aS oy -E . : : 
of this statistic are approximated and the distribution of Sa Ela} is approximated with 


¢ 
& 
a Student’s t distribution. The desired confidence interval is constructed using the ap- 


proximated Student’s t distribution. The analysis in the following paragraph serves only 


as a means te find a plausible expression for the degrees of freedom of the approximated 


Student's distribution. 








dt 





x-—7 
oe denotes a value constructed from 
lsh 5} 


observed data. The upper case version of the same expression denotes a random Vvari- 


We begin as in equation (4.4), where 


be able. 


ie 


l-a=P 


al ofn + o?, ne Gens xX-p af stin + sofm 
/ > ~ $4 f > 
nN Ss? + Sy af Sy]n + st] I st + sy 


—== 
af oxln + oy lm 


X — ¥ = (uy — wy)-+ (ey — ay) 





Ly — By 
| Ze yx al 
» 4 
foyla + oyjm x-y s 
=p i aN ed nae < ae Le Se , (4.5) 
4 2 2 
S$ + S$ syin + splm aj sxln + s $]m 
otln + ohn sh + st 


where o3, 3, s3, and s? denote constants. The last probability in Equacion (4.5) suggests 


— 


-4 noncentral t distribution. Consequently the random variable 
[ (sin + shmys? + 5} 
V (cin + opm) st + s?) 
in the denominator should be a random variable of the form y ae That is, 


2 
sdjn + s}fm Ss} + Sy 2 
SS OS = 


2 


ohn + oF] s+ s} 





. : ee 264, ‘ 2G}, 
Since Var(y2) = 2v , Var(S}) = ae a and Var(S}) = Poa Be 
. 2 stin + sim 2 ] 26% 2e}: 
Re er oe EY en) 
SH S$ (oy[n + oyfmye \ PO mi 


Thus, 














hash \, (ohn + ohm)? 
yo ap aceite Teguee Ad ST eA A 
s3]n + ffm ; of] (n N+ oyl(n —1) 


2 ( Stest \p__ (ohne shy 
: sdjn + s3fm sla -l)+ spin -1) 


2, 22 
(sy + Sy) 


Syl(n —1) + sy/(m — 1) 


In summary, if we fit a t distribution to the distribution of ~Ele) , the expression in 
- p 
3 


the Equation-(4.6) might be a plausible random value for the degrees of freedom. The 


results of computer simulations: will indicate the accuracy of this.approximation. 


In order to formulate the lower confidence interval statement.for = , Wwe first 


ach + oF 
. We then ap- 


A= 
[Sh 5 


~“N 
with a central Student’s t distribution with v-degrees 


find the mean, E{g]=4y,, and een Var[g] = 03, for g= 





proaimate the distribution of 2 
of freedom, where v is given in sacaden (4.6). The lower confidence limit, 2, ya, £08 


u, Will yield the corresponding lower confidence limit for reliability 
i A 
Ry, jaa = Oli, ,_,) 


We now proceed to find yw, and o?. The Taylor expansion of g at (uy — Hy, o} + 4) , Us- 


ing only first order derivatives is given by: 


g( — ¥, Se+ S}) = slay — ayo +09) + LY — ¥ - (ux — ay) 





(Ly Hy, of io o}) 





ae aS 
+ (5% 453-(0 + +o) 2 2+R 3 
Fs x Y. as? +S. 3) (Hy - Hy, of 04) fa (4. ) 
mm py X—-¥ (uy 3 — 4 
ve EL. ee lux Hy) —Es3+ S}- (c4-+ o})] 2h + R, 
Jere — foh+eh eh + o4y" 


where &, includes terms that converge to 0 at the same rate as max( + ; t ie as n and 
m become large. 


The expected value and variance of g are 
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a EXT HY 
of + of 
and 
in Sue 
Var[g] = 02 
I =, Wyn my)? 2,2 
= ar Sears Var(¥ ~— ¥) + Se ae Var(oy + o}) 
Cyt cy (oy + oy)” 
4 4 
ee otfn --o%]m (Hy —py) ( 2¢y 2c} ) 
ort ot Aoi +cz) \2al° m-i 
_ ant ofl 4 lew J ok int +) 
Sloped Norte g2ty \ melo mal | 
Cyt oy (oy + oy) 


An estimator for o, is 


» _| Stine Shim | (¥-¥P (Se SFY ip 
os= Foes 2, ¢2y3 Fol ae : 
Sy t Sp 2(S¥ + S¥) 


Then an approximate confidence limit for y, is the follows: 


A =e A 
Mey, lea g @ ghi_n, ¥v 


A 
= OL, lax 


(4.8). 


(4.9) 


(4.10) 


(4.11) 


The corresponding 100(1 — «) % lower confidence limit, Riga, , for component reliability 


is 
Ry, I-x = (8), ma) 


1. Example 


(4.12) 


Let x, denote the strength of solid missile motor chambers that were pressurized 


until they burst, i= 1, 2, .., 25. Let 3, denote the maximum pressure observed on 30 


solid missile motors that use this type of chamber. X and Y are assumed to have normal 


distributions with unequal variances. The coded X and Y data are: 


25 














X: 292.65 301.86 297.86 312.67 300.61 292.18 297.55 


284.64 297.93 295.22 315..49 294.10 299.85 292. 24 
312.25 298.87 295.41 292.41 287.71 298.57 300. 34 


303.28 309.52 283.32 308.52 


¥: 242.52 261.00 229.37 261.:82 248.36 174.04 192.50 
265.64 199.70 225.52 237.43 285.66 199.20 203.65 
274.93 302.13 302.85 233.59 275.95 220.46 318.34 
269.86 200.10 318.19 216.69 281.78 301.20 348.56 


234.15 234.45 


The corresponding 90% lower confidence li:nit for component reliability, 


R= P(X > ¥) is computed’ as follows: 
n=25, ¥= 298.60, s#= 70.25 
m= 30, p=251.99, s%= 1873.79 


y.,di9 
(sy + 55) 


ar a a = 3) 
Sy](t — 1) + sy](m — 1) 


2 a 4 

A syln + sy]m (R= 7P SY sy 1 

Cee ee ee aN Bool tea i) ee 
S¥ + SY 2(s¥ + Sy) u 7 


e A, - 
51,0.9 = & — &sfo,9, 31 = 0.76 


Rios = (6, 0.9) = 0.78 

















The corresponding 90%. Govindarajulu nonparametric confidence limit is 0.69. This is 
not surprising because our procedure uses more information about the distributions of 
& and Y. Intervals .stimated by nonparametric procedures are usually wider than those 
estimated by parametric procedures. The amount of the difference is somewhat sur- 


prising-for these sample sizes of 25 and 30. 


B. ACCURACY OF THE PROCEDURE 
1. Computer simulation 
a. Simulatien procedure 
Computer s:mulations were used to determine its accuracy for specific sets 
values of 1, nl, Uy, Hy, 62, and 3}. 
Parameters are chosen in a Way to cover practical conditions. The sets of 


parameters are as follows: 
© Uy> By 
¢ Reliability, R : 0.90, 0.95, 0.99 
¢ Standard Deviations 


Oy: TQ 10.0 
oy: 2:0, 40.0 


¢ Sample sizes, (n,m): (10, 20); (25, 35); (75, 50) 
* Confidence leyzi, 2 + 0.05, 0-10, 0.20 . 


Actual sets of parameters are tabulated in Appendix E. 





Simulation procedure is summarized as follows: 
p 


1) Generate n random norimal variates X, 1 ~ \(u,6%):; m random normal variates 
Y, ¥~ N(uy63). 


2) Compute J, ¥,, $2. 5}. 
3) Compute 6,6, ;-,. 


4) Compute £1, = 0(6,,1.,), for = 0.05, 0.10, 0.20. 


5) Repeat steps | through 4, for 1000 times. Then order the R,)., to get 
Ray taesy Rentgeay ot Riesa-sta-sy 9» Race. ta-2» from the smallest to the largest. 











6) Print Rios), 20-3): 
7) Print max{Ro,ra-ay' Ro,za- < RX} the actual confidence level of the -approximation 


will-be 1000" 


8) Compute and print the sample mean.and variance-of the estimation error. 


b. Simulation language 
The simulations were conducted on the N.P:S. mainframe IBM 3033 com- 
puter. The-programming language used for the simulations is VS FORTRAN 2. The 
random number generator LLRANII was used to generate normal variates. The IMSL 
Statistics function TIN was used to compute the percentile of Student’s.t distribution, 
and ANORDF was used to compute the probability of standard normal distribution. 
The FORTRAN source-code for the simulation of the developed approximate procedure 
are attached in Appendix D. The FORTRAN code for simultaneous -comparison be- 
tween the approximate procedure and nonparametric procedure are attached in Appen- 
dix F. 
2. Analysis of simulation results 
The simulation results are tabulated in Tables § through 10. The results show 
the developed approximate procedure is quite accurate, because for every case the 
100(1 — «)" percentile point of R., are all very close to the true reliability. Also, the 
‘mean error from R’ and ‘Variance of error’ reduce rapidly with the increases of sample 
size. Comparison simulations were run for the approximate procedure versus the 


Govindarajulu nonparametric procedure. Results are tabulated in table 1] and 12. 
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Table 8. as hepa OF 80% CONFIDENCE LIMIT APPROXIMATION OF 
RE PIA> Y) FOR UNEQUAL VARIANCES CASE 


True confi- | Meanerror | Variance of | 
Root», t-2) | dence level from R error 
=a 3008 
20 3996 021 | 
, Psat er [sr | 0011 
~7310 0082 
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Table 9. ANALYSIS OF 90% CONFIDENCE LIMIT APPROXIMATION OF 
R = P(X> Y) FOR UNEQUAL VARIANCES CASE : 
Mean error | Variance of 


: ee * True-confi- 
x } Ricoot-v, -2) | dence level from R error 


9008 8970 0893 
2.0 
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‘0013 
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0032 
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0132 0001 
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Table 10: pie OF 95% CONFIDENCE LIMIT APPROXIMATION OF 
3 R= P(X> ¥) FOR UNEQUAL VARIANCES CASE 
Variance of 


: True-confi- | Mean error 
Rani ux | dence level from R error 
9560 A218 | .0059 


25 ‘8988 9530 0816 0027 


75| 0 8998 
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Table 11. COMPARISON OF DIFFERENT PROCEDURES ON 90% LOWER 
CONFIDENCE LIMIT ESTIMATION OF R = P(X> Y), UNEQUAL 
eee CASE | 


jor) er || : Approximate Estimation _Nonparametri Bound — |- 
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10001-7). L(i=2) 100%t~2), L(i=7) 
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Table 12. COMPARISON OF DIFFERENT PROCEDURES ON 95% LOWER 
CONFIDENCE LIMIT ESTIMATION OF R = P(X> Y), UNEQUAL 
; VARIANCES: CASE 
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V. CONCLUSIONS. AND RECOMMENDATIONS 


The-lower confidence limit estimation procedures developed-in this thesis for equal 
variances case as well-as for unequal variances cases, are very accurate. These proce- 
dures are simple to evaluate and require only the use of central- Student’s t tables, in 
contrast to the existing parametric procedures of this type which require extensive use 
of noncentral Student’s t tables. 


Although these procedures are developed for lower confidence limits, upper-or two- 


sided confidence limits -for the reliability are readily obtained. 








; APPENDIX A. FORTRAN CODE FOR INTERVAL ESTIMATION 
PROCEDURE - NORMAL EQUAL VARIA NCES -CASE 


, PROGRAM EQSIGM 
seavdeieikdedeiedesetcascdedsdeseacdedeveredsdeedcdevcdciederevedeadesededededeoderedederedscaedede de aese dea sek sede ve de ag se vere Fee 
ve a 
* THIS PROGRAM IS TO VALIDATE THE LOWER CONFIDENCE BOUND APPROXI- * 
* MATION PROCEDURE FOR P( X > Y ),.-WHERE X, Y ARE NORMALLY DIS- * 
*  TRIBUTED WITH UNKNOWN MEANS AND A UNKNOWN BUT EQUAL VARIANCE. % 
* ve 
* VARIABLES DESCRIPTION: s 
. ALPHA - NOMINAL CONFIDENCE LEVEL * 
¢ ANORDF - IMSL FUNCTION FOR NORMAL PROBABILITY * 
7 CASE - NUMBER OF TEST PARAMETER SETS % 
re ER - ERROR BETWEEN LIMIT AND TRUE RELIABILITY * 
* ERBAR - AVERAGE OF ER : * 
¢ ERSSQ - SUM OF SQUARES OF ER * 
° ERSUM - §UM OF ER * 
‘ ERSV - SAMPLE VARIANCE OF ER * 
? CLOSE - INDEX OF THE CLOSEST ESTIMATE : * 
% DELTA - NONCENTRALITY OF T DISTRIBUTION * 
* DF - DEGREE OF FREEDOM OF T DISTRIBUTION * 

‘ * K ~ STATISTIC TO ESTINATE DELTA * 
* LNORM - RANDOM NUMBER GENERATOR FOR NORMAL VARIATES * 
¢ M - SAMPLE SIZE OF Y RANDOM VARIABLE % 

* ? MUX - POPULATION MEAN OF X RANDOM VARIABLE * 
¥ MUY - POPULATION MEAN OF Y RANDOM VARIABLE ¥ 
? N - SAMPLE SIZE OF X RANDOM VARIABLE * 
* R - REAL RELIABILITY * 
* REP - REPETITION OF SIMULATION * 
te RLHAT - LOWER CONFIDENCE LIMIT OF RELIABILITY * 
* SIGMAX - POPULATION STANDARD DEVIATION OF X * 
* SIGMAY - POPULATION STANDARD DEVIATION OF Y % 
* SP - POOLED SAMPLE VARIANCE OF X AND Y a 
ve SUMSQX - SUM OF SQUARES OF X * 
w SUMSQY - SUM OF SQUARES OF Y * 
* SUMX - §SUM OF X % 
” SUMY - SUM OF Y i 
* TIN - IMSL FUNCTION TO COMPUTE PERCENTILE OF T DIST. % 
* XBAR ~ AVERAGE OF X * 
ve YBAR - AVERAGE OF Y * 
sy Se 
% REQUIRED EXTERNAL FUNCTIONS: LNORM, TIN, ANORDF * 


Jeni atvalpatvaleaieatenipeion anincteuteaivniaatenten'pateatulanton wateateateateuma's aia seal 


ale 
7 


INTEGER REP, CASE 
REAL ALPHA, TWO 
PARAMETER( ALPHA=0. 05) 
PARAMETER( TWO=2. 0) 


w 
A 











” 


200 


ve 


PARANETER( REP=1000) 
PARAMETER( CASE=18) 


INTEGER I, 11, J, U, V, XSEED, YSEED, N, M, CLOSE 

REAL MUX, MUY, XBAR, YBAR, SIGMAX, SIGMAY, DF, R, RLHAT(REP), 
+ K, DELTA, XX(100), YY(100), X, Y, SUMSQX, SUMSQY,. SUMX, 
s SUMY, TIN, ANORDF, RN, RM, DIFF, SP, TEMP, ER, ERSUM, 
+ ERSSQ, ERBAR, ERSV 


CALL EXCMS('FILEDEF 12 DISK SETUP1 DATA Al') 
CALL EXCMS('FILEDEF 18 DISK OPT1 DATA Al’) 


DO 2000 I=1, CASE 
READ (12,2200) XSEED, YSEED, MUX, MUY, SIGMAX, SIGMAY, N, M, R 
DF = REAL(N+M-2) 
RN = REAL(N) 
RM = REAL(M) 
ERSUN = 0.0 
ERSSQ = 0.0 
DO 1000 J=1, REP 
CALL LNORM(XSEED, XX, N, 2, 0) 
CALL LNORM(YSEED, YY, M, 2, 0) 
SUNSQX = 0.0 
SUMSOY = 0.0 
SUMX = 0.0 
SUMY = 0.0 


< TRANSFORM X, Y TO DESIRED PROPERTIES > 


DO 200 U= 1, N 
= XX(U) * SIGMAX + MUX 


SUNSQX= SUMSQX + X * X - 
SUMX = SUNX + X 
CONTINUE 


XBAR= SUMX / RN 

DO 300 V=1, M 
Y = YY(V) * SIGMAY + MUY 
SUMSQY = SUMSQY + Y * Y 
SUMY = SUMY + Y 

CONTINUE 


< COMPUTE CONFIDENCE LIMIT FOR RELIABILITY > 


YBAR = SUMY / RM 

SP = SQRT( (SUMSQX - RN“XBAR*XBAR + SUMSQY - RM*YBAR*YBAR) 
+ / DF) 

K = MAX( (XBAR - YBAR) / SP, 0.0 ) 

DELTA = K - SQRT( (RN+RM)/(RN*RM) + K*K / (2. 0%*(RN+RM-2. 0))) 
rs ** TIN(1. 0-ALPHA, DF) 

RLHAT(J) = ANORDF(DELTA/SQRT(TWO)) 


< COMPUTE THE MEAN AND VARIANCE OF ER > 


ER = R - RLHAT(J) 
ERSUM = ERSUM + ER 











ERSSQ = ERSSQ + ER * ER 
1000 CONTINUE 
ERBAR = ERSUM / REALCREP) 
ERSV = ( ERSSQ - REAL(REP) * ERBAR * ERBAR ) / REALCREP-1). 


¥ -< SORT CONFIDENCE LIMITS IN ASCENDING ORDER > 


DIFF = 2.0 
DO 1800 Ti=1, REP 
DO 1500 J=Il+1, REP 
IF (RLHAT(J) .LT. RLHAT(I1)) THEN 
TEMP = RLHAT(11) 
RLHAT(I1) = RLHAT(J) 
-RLHAT(J) = TEMP 
ENDIF 
1500 CONTINUE 


* < FIND THE CLOSEST CONFIDENCE LIMIT ESTINATE > 


IF ((CR-RLHAT(I1)) .GE. 0.1E-6) . AND. 
+  ((R-RLHAT(I1)) .LE. DIFF)) THEN 
DIFF = R - RLHAT(I1) 
CLOSE = Il 
_ ENDIF 
1800 CONTINUE 
RITE (18,2100) I, MUX, N, SIGMAX, MUY, M, SIGMAY, R, 
RLHAT(NINT(REAL(REP)**(1.0-ALPHA))), RLHAT(CLOSE), 
REAL(CLOSE)/1000..0, ERBAR, ERSV 
2000 “coNTINUE 


2100 FORMAT(' SIMULATION: ',12,/,'MUX: ' sF5.1,T16, We *ot2y 
T35,'SIGMAX: ' FA. 1 »/,\ MUY: PIS; T16, 'm: ',12, 

T35,'SIGMAY: '.F4. 1 s/) 'TRUE R': ',F7. 5, 

T35,'RLHAT: ',F7.5,/,' CLOSEST RLHAT: ',F7.5, 

735, ‘TRUE CONFIDENCE LEVEL: | ,F5.3,/, 

‘MEAN ERROR WIDTH: ' ,F7.5,/, 

"VARIANCE OF ERROR: ' ,F7.5,///) 

2200 FORMAT (I5,1X,15,1X,F5.1,1X,F7.3,1X,F4. 1,1X,F4.1,1X,12,1X,12, 1X, 


bbb bb + 























APPENDIX B. 


93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 
93943 


COCO OCoO OOOO OOCCCOOCDSD 


x 
MEAN 


SIMULATION PARAMETER SETS FOR EQUAL 


8. 187 
8. 187 
8. 187 
63. 740 
63. 740 
63. 740 
7.674 
7,674 
7.676 
53.472 
53.472 
53.472 
6.711 
6.711 


Y 
MEAN 


VARIANCES CASE 


1.0 1.0 
1.0 1.0 
1.0 1.0 
20.0 20.0 
20.0 20.0 
20.0 20.0 
1.0 1.0 
1.0 1.0 
1.0 1.0 
20.0 20.0 
20..0 20.0 
20. 0 20.0 
1.0 1.0 
1.0 1.0 
1.0 1.0 
20.0 20.0 
20.0 20.0 
20.0 20.0 


X ¥ 
STDV STDV 


8 
8 


20 


8 
30 
30 

8 
30 
30 

8 
30 
30 

8 
30 
30 


: SETUP] DATA (FOR APPROXIMATE PROCEDURE) 


. 900 
900 
900 
900 
900 
. 900 
950 
950 
950 
950 
950 
950 
990 
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FILE: SEQ4 DATA (FOR COMPARISON SIMULATIONS) 


16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 


COSDDOCDDDDCDDOCOSOONSO 


SEED SEDD 


MEAN 


. 187 
. 187 
. 187 
63. 740 
63. 740 
63. 740 
7.674 
7.674 
7.674 
53.472 
53.472 
53.472 
6.711 
6.711 


ooo 


ID hb 


PER OOORPEHOOO rH PE 


0 


Nh bo 
ooooococooooooceoo 
to 
9 OE ES Sea ees 
oooooooo°o°o°ededo°o°o°ecoo 


25 
50 
90 


25- 
50 


MEAN 


STDV STDV 
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APPENDIX -C. FORTRAN CODE FOR SIMULTANEOUS COMPARISON 
SIMULATION OF APPROXIMATE PROCEDURE VS. 


NONPARAMETRIC PROCEDURE - EQUAL VARIANCES 
PROGRAM COMEQ 


devededevevesedesetedededevervesededevedesese teveve terete retest Fe Ke Fe Ke a0 FTE Tee TENET TENET TTL T ER NCTE TNC TCT TON NET TEETER TC INE 


* THIS PROGRAM IS TO COMPUTE THE LOWER CONFIDENCE LIMIT OF % 
* R= P( X > Y ), WHERE X, Y ARE INDEPENDENTLY NORMALLY DIS- * 
* TRIBUTED WITH UNKNOWN MEANS AND A UNKNOWN BUT EQUAL VARIAN- * 
* CES WITH APPROXIMATE PROCEDURE AND WITH NONPARAMETRIC PRO- * 
¥ CEDURE * 
wv * 
* VARIABLES DESCRIPTION: * 
1 * 
a -ALPHA - NOMINAL CONFIDENCE LEVEL * 
* ANORDF - IMSL FUNCTION FOR NORMAL PROBABILITY * 
* ANORIN - IMSL FUNCTION FOR INVERSE NORMAL CDF * 
* CASE - NUMBER OF TEST PARAMETER SETS * 
% BIGU - MANN-WHITNEY STATISTIC bd 
% CLOSE ~ INDEX OF THE CLOSEST ESTIMAT (APPROXI PROCEDURE) * 
* CLOSEN - INDEX OF THE CLOSEST ESTIMAT (NONPARA PROCEDURE) * 
¥ DELTA - NONCENTRALITY OF T DISTRIBUTION ze : 
bd DF - DEGREE OF FREEDON OF T DISTRIBUTION * 
* EPS ~ WIDTH OF THE CONFIDENCE BOUND * 
% EROR - EROR OF ESTIMATION (APPROXI PROCEDURE) * 
% ERORN - EROR OF ESTIMATION (NONPARA PROCEDURE) ¥ : 
* ERBAR - MEAN OF EROR (APPROXI PROCEDURE) # 
¥ ERBARN - MEAN OF EROR (NONPARA PROCEDURE) as 
ve ERSSQ ~ SUM OF SQUARE OF EROR (APPROXI PROCEDURE) A 
ve ERSSQN - SUM OF SQUARE OF EROR (NONPARA PROCEDURE) * 
ERSUM - SUM OF EROR (APPROXI PROCEDURE) * 
= ERSUMN - SUM OF EROR (NONPARA PROCEDURE) * 
v ERSV ~ SAMPLE VARIANCE OF EROR (APPROXI PROCEDURE) * 
a ERSVN ~ SAMPLE VARIANCE OF EROR (NONPARA PROCEDURE) * 
ve K - STATISTIC TO ESTIMATE DELTA % 
ae LNORM - RANDOM NUMBER GENERATOR FOR NORMAL VARIATES % 
a M - SAMPLE SIZE OF Y RANDOM VARIABLE * 
# MUX - POPULATION MEAN OF X RANDOM VARIABLE * 
* MUY - POPULATION MEAN OF Y RANDOM VARIABLE = 
ae N - SAMPLE SIZE OF X RANDOM VARIABLE * 
ve NU ~ THE SMALLER OF SAMPLE. SIZES * 
oo R ~ REAL RELIABILITY * 
ve RB - CONFIDENCE BOUND OF THE RELIABILITY * 
ae REP - REPETITION OF SIMULATION = 
a RLHAT ~- LOWER CONFIDENCE LIMIT OF RELIABILITY * . 
% RTILD ~ POINT ESTIMATOR OF THE RELIABILITY = 
¥ SIGNAX - POPULATION STANDARD DEVIATION OF X * 
* SIGMAY - POPULATION STANDARD DEVIATION OF Y = 
bad SP - POOLED SAMPLE VARIANCE OF X AND Y = : 
= SUMSQX - SUM OF SQUARES OF X cs 
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200 
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SUMSQY - SUM OF SQUARES OF Y * 


SUMX - SUM OF X ve 

SUMY - SUM OF Y * 

TIN - IMSL FUNCTION TO COMPUTE PERCENTILE OF T DIST. * 

XBAR - AVERAGE OF X * 

YBAR - AVERAGE OF Y * 

ve 

REQUIRED EXTERNAL FUNCTIONS: LNORM, TIN, ANORDF, ANORIN * 


INTEGER REP, CASE 
REAL ALPHA, TWO | 
PARAMETER( ALPHA=0. 05) 
PARAMETER(TWO=2..0) 
PARAMETER( REP=1000) 
PARANETER( CASE=18) 
INTEGER I, Ii, J, U, V, XSEED, YSEED, N, M, CLOSE, 
A, B, Al, Bl, CLOSEN 
REAL MUX, MUY, XBAR, YBAR, SIGMAX, SIGMAY, DF, R, RLHAT(REP), 
é K, DELTA, X(120), Y(120), X1(120), Yi(120), SUMSQX, 
+ SUMSQY, SUMX, SUMY, TIN, ANORDF, RN, RM, DIFF, SIGMA, TEMP, 
+ ER, ERSUN, ERSSQ, ERBAR, ERSV, 
+ BIGU, RTILD, NU, EPS, ANORIN,ERN, ERSUMN, ERSSQN, ERBARN, 
+ ERSVN, DEF, RB(REP) 


Se 


ee ee ee ee eee eee) 


CALL EXCNS('FILEDEF 12 DISK SEQ4 DATA A1') 
CALL EXCNS('FILEDEF 18 DISK AEQ4 DATA Al') 
DO 2000 I=1, CASE 
READ (12,100) XSEED, YSEED, MUX, MUY, SIGMAX, SIGMAY, N, M, R 
FORMAT (15,1X,15,1X,F5. 1,1X,F7.3,1X,F4. 1,1X,F4. 1,1X,12,1X,12,1%, 
+5. 
DF = REAL(N+N-2) 
RN = REAL(N) 
RM = REAL(M 
ERSUNN = 0.0 
ERSSQN = 0.0 
ERSUN = 0.0 
ERSSQ = 0.0 
DO 1000 J=1, REP 
CALL LNORM(XSEED, X, N, 2, 0) 
CALL LNORM(YSEED, Y, M, 2, 0) 
SUMSQX = 0.0 
SUMSOY = 0.0 
SUMX = 0.0 
SUNY = 0.0 
DO 200 U= 1, N 
X1(U) = X(U) * SIGMAX + MUX 
SUNSOX= SUMSQX + X1(U) * X1(U) 
SUMX = SUNX + X1(U) 
CONTINUE 
XBAR= SUMX / RN 
DO 300 V=1, M 
Y1(V) = Y(V) * SIGMAY + MUY 
SUMSQY = SUMSQY + Y1(V) * Y1(V) 


Vou tv 
—_ 
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SUMY = SUMY + Y1(V) 
300 CONTINUE 


* PROCEDURE FOR PARAMETRIC 

YBAR = SUMY / RMN 

SIGMA. = SQRT( (SUMSQX - RN**XBAR*XBAR + SUMSQY - RM*YBAR*YBAR) 
+ _ { DF) : 

K = MAX( (XBAR - YBAR) / SIGMA, 0.0 ) 

DELTA = K - SQRT( (RN+RM)/CRN*RM) + K*K / (2. O*(RNtRN-2..0))) 
+ * TIN(1. 0-ALPHA, DF) 

RLHAT(J) = ANORDF(DELTA/SQRT( TWO) ) 

ER = R - RLHAT(J) 

ERSUM = ERSUN + ER 

ERSSQ = ERSSQ + ER * ER 


* PROCEDURE FOR NONPARAMETRIC 

BIGU = 0.0 

DO 500 A=1, N 

DO 400 B = 1, M 
IF (X1(A) .GT. Y1(B)) BIGU = BIGU + 1.0 

400 CONTINUE 
500 CONTINUE 

RTILD = BIGU / (RN * RM) 

NU = MIN (RN, RM) 

EPS = 1.0 / SQRT(4.0 * NU) * ANORIN(1.0 - ALPHA) 

RB(J) = RTILD - EPS 

ERN = R - RB(J) 

ERSUNN = ERSUMN + ERN 

ERSSQN = ERSSQN + ERN * ERN 
1000 CONTINUE 


ERBARN = ERSUNN / REAL(REP) 
ERSVN = ( ERSSQN ~- REAL(REP) * ERBARN * ERBARN ) / REAL(REP-1) 
DFF = 2.0 
DO 1300 Al = 1, REP 
DO 1200 Bi = Ail + 1, REP 
IF (RB(B1) . LT. RB(Ai)) THEN 
TEMP = RB(A1) 
RB(A1) = RB(B1) 
RB(B1L) = TEMP 
ENDIF 
1200 CONTINUE 
IF ((CR-RB(A1L)) .GE. 0.1E-6) . AND. 
+  ((R-RBCAl1)) .LE. DFF)) THEN 
DFF = R - RBCA1) 
CLOSEN = Al 
ENDIF 
1300 CONTINUE 
Naw ee EM EO RE NE EW 
* PROCEDURE OF PARAMETRIC 
ERBAR = ERSUN / REALCREP) 
ERSV = ( ERSSQ - REAL(REP) * ERBAR * ERBAR ) / REAL(REP-1) 
DIFF = 2.0 
DO 1800 Ii=1, REP 
DO 1500 J=Ii+1, REP 
IF (RLHAT(J) . LT. RLHAT(I1)) THEN 















1500 


1800 


1900 


1910 
1920 


2000 





+ T16, ‘SIGMA : F4, 1) 


+ F5.4,//) 


TEMP = RLHAT(I1) 
RLHAT(11) = RLHAT(J) 
_ RLHAT(J) = TEMP 
ENDIF 
CONTINUE 
IF (((R-RLHAT(I1)) .GE. 0.15-6) . AND. 
((R-RLHAT(I1)) .LE. DIFF)) THEN 
DIFF = R - RLHAT(I1) 
CLOSE = I1 
ENDIF 
CONTINUE 
WRITE (18,1900) I, N, M, R, SIGMAX 
WRITE (18,1910) RLHAT(NINT(REAL(REP)**(1. -0-ALPHA))), ERBAR 
WRITE (18,1920) RB(NINTCREAL(REP)*( 1. 0- mabe es a 
FORMAT(’ CASE: ,12 o/s! N: ', 12, T16, 'M: ', 12, /, 'R: ', F4.3, 


FORMAT('< PARAMETRIC ar /; * RLHAT: ', F5.4, T16, ‘ERROR: ', F5.4) 
FORMAT('< NONPARAMETRIC >', /, 'RLHAT: ', F5.4, T16, ‘ERROR: ', 


CONTINUE 
STOP 
END 




















APPENDIX D. FORTRAN CODE-FOR INTERVAL ESTIMATION 


PROCEDURE - NORMAL UNEQUAL VAR IANCES CASE 


PROGRAM UNEQSN 
* 
Tetedede deere dees Te TTI I Le TCN TE HL TI TE TE TENT IE TE TE ETE TE TN LETTE TE ENE TENE NEEDED TE TELE TE TENE NTE HEE NE TE TE TET TE TE 
* THIS PROGRAN IS TO VALIDATE THE LOWER CONFIDENCE LIMIT APPROXI- * 
* MATION PROCEDURE. FOR P( X > Y )-, WHERE X, Y ARE INDEPENDENTLY * 
* NORNALLY DISTRIBUTED WITH UNKNOWN MEANS AND UNKNOWN AND UNEQUAL * 


%* — VARTANCES * 
cy oe 
% ALPHA - NOMINAL CONFIDENCE LEVEL ¥ 
% ANORDF - IMSL FUNCTION FOR NORMAL PROBABILITY * 
ve ASVX - SAMPLE VARIANCE OF X DEVIDED BY SAMPLE SIZE ¥ 
* ASVY - SAMPLE VARIANCE OF Y DEVIDED BY SAMPLE SIZE % 
* CASE - NUMBER OF TESTING ve 
¥ ERROR - ERROR OF ESTIMATION 
ag ERRBAR - MEAN OF ERROR ¥ 
¥ ERRSSQ - SUM OF SQUARE OF ERROR * 
ve ERRSUM - SUM OF ERROR ve 
¥ ERRSV - SAMPLE VARIANCE OF ERROR 7 
3 CLOSE - INDEX OF THE CLOSEST ESTIMAT * 
: DELTAH - DELTAHAT; A ESTIMATOR OF DELTA 3 
¢ DF - DEGREES OF FREEDON OF DELTAHAT % 
LNORM - RANDOM NUMBER GENERATOR FOR NORMAL VARIATES % 
3 N - SAMPLE SIZE OF RANDOM VARIABLE X % 

MUX ~ POPULATION MEAN OF X a 
€ MUY - POPULATION MEAN OF Y * 
? N - SAMPLE SIZE OF RANDOM VARIABLE Y ve 
° R - REAL RELIABILITY ve 

REP - REPETITION OF SIMULATIONS ¥ 
* RLHAT - LOWER CONFIDENCE LIMIT OF RELIABILITY ve 
* SDHAT ~ SAMPLE STANDARD DEVIATION OF THE DELTAHAT * 
* SIGMAX ~ STANDARD DEVIATION OF X ¥ 
* SIGHAY - STANDARD DEVIATION OF Y a 

SUNSQX - SUM OF SQUARE OF RANDOM SAMPLE OF X 7 

SUNSQY - SUM OF SQUARE OF RANDOM SAMPLE OF Y * 
* SUMX - SUM OF RANDOM SAMPLE OF X ve 
* SUNY - SUN OF RANDOM SAMPLE OF Y 
* SVX - SAMPLE VARIANCE OF X * 
* SVY - SAMPLE VARIANCE OF Y isd 
* T - PERCENTILE OF THE T DISTRIBUTION ve 
ve TIN - IMSL FUNCTION TO COMPUTE PERCENTILE OF T DISTRIBUTION * 
ve VARHAT ~- SAMPLE VARIANCE OF DELTAHAT ve 
XBAR - SAMPLE MEAN OF X * 
* YBAR - SAMPLE MEAN OF Y ve 


% REQUIRED EXTERNAL FUNCTIONS: LNORM, TIN, ANORDF = 


ale 
o 





Se we om oo oe 


100 


200 





INTEGER REP, CASE 

REAL ALPHA 

PARAMETER( REP=1000) 

PARAMETER( CASE=18) 

PARAMETER( ALPHA=0. 05) 

INTEGER I, T1, J, U, V, XSEED, YSEED, N, M, CLOSE 

REAL RM, RN, MUX, MUY, XBAR, YBAR, SIGMAX, SIGMAY, R, 
+ X(100), Y¥(100), X1, Y1, SUMSQX, SUMSQY, SUMX, SUMY, 
+ SVX, SVY, ASVX, ASVY, DF1, DF, T, TIN, DELTAH, VARHAT, 
+ SDHAT, RLHAT(REP), ANORDF, TEMP, DIFF, ERROR, ERRSUM, ERRSSQ, 
+ ERRBAR, ERRSV 

CALL EXCMS('FILEDEF 12 DISK SETUNQ DATA A1') 

CALL EXCMS('FILEDEF 18 DISK OPT2 DATA Al‘) 

WRITE (18,%) ‘UNEQUAL VARIANCES’ 

DO 2000 I=1, CASE 

READ (12,100) XSEED, YSEED, MUX, MUY, SIGMAX, SIGMAY, N, M, R 
FORMAT (15,1X,15,1X,F5. 1,1X,F7. 3,1X,F4. 1,1X,F4. 1,1X,12,1X,12,1X, 
+F5. 3) 

RN. = REAL(N) 

RM = REAL(M) 

ERRSUMN = 0. 0 

ERRSSQ = 0.0 





< TRANSFORMATION OF X, Y TO DESIRED PROPERTIES > 


DO 1000 J=1, REP 

CALL LNORM(XSEED, X, N, 2, 0) 

CALL LNORM(YSEED, Y, M, 2, 0) 

SUMSQX = 0.0 

SUNSOY = 0.0 

SUMX = 0.0 

SUNY = 0.0 

DO 200 U= 1, N 
X1 = X(U) * SIGMAX + MUX 
SUMSQX= SUMSQX + X1 * Xi 
SUNX = SUMX # XI 

CONTINUE 

XBAR= SUNX / RN 

DO 300 V=1, M 
Y1 = ¥(V) * SIGMAY + MUY 
SUNSQY = SUMSQY + Y1 * Y1 
SUMY = SUMY + Yi 

CONTINUE 


< COMPUTE CONFIDENCE LIMIT OF RELIABILITY > 
YBAR = SUMY / Ri 

SVX = (SUMSQX - RN * XBAR * XBAR) / (RN - 1.0) 
SVY = (SUMSQY - RN * YBAR * YBAR) / (RM - 1.0) 


ASVX = SVX / RN 

ASVY = SVY / RY 

Fl = (SVX + SVY) * (SVX + SVY) / ( SVX*SVX / (RN-1.0) + 
+ SVY*SVY / (RM-1.0) ) 


DF = ANINT(DF1) 
T = TIN( 1.0-ALPHA, DF) 
DELTAH = (XBAR - YBAR} / SQRI(SVX + SVY) 














VARHAT = (ASVX + ASVY) / (SVX + SVY) + 


+ (XBAR - YBAR) * (XBAR - YBAR) / 
+ ( 2.0*(SVX + SVY)%*3 ) % 
+ ( -CSVX**2) / ERN-1.0) + (SVY%"2)- / (RM-1.0) ) 
SDHAT = SQRT(VARHAT) 
RLHAT(J) = ANORDF(DELTAH - T * SDHAT) 
* 
*: < COMPUTE THE MEAN AND VARIANCE OF ESTIMA’. ‘N ERROR > 


ERROR= R ~ RLHAT(J) 


ERRSUN = ERRSUM + ERROR 
ERRSSQ = ERRSSQ + ERROR * ERROR 
ERRBAR = ERRSUM / REALCREP) 


ERRSV = ( ERRSSQ - REAL(REP) * ERRBAR * ERRBAR ) / REAL(REP-1) 
1000 CONTINUE 


* < SORT CONFIDENCE LIMITS IN ASCENDING ORDER > 


DIFF = 2.0 
DO 1800 Ii=1, REP 
DO 1500 J=I1+1, REP 
IF (RLHAT(J) .LT. RLHAT(I1)) THEN 
TEMP = RLHAT(11) 
RLHAT(1I1) = RLHAT(J) 
RLHAT(J) = TEMP 
ENDIF 
1500 CONTINUE 


a < FIND THE CLOSEST CONFIDENCE LIMIT ESTIMATE > 


IF (((R-RLHAT(1I1)) .GE. 0.1E-6) . AND. 
+  ((R*RLHAT(I1)) .LE. DIFF)) THEN 
DIFF = R ~ RLHAT(I1) 
CLOSE = Il 
ENDIF 
1800 CONTINUE 
WRITE (18,1900) I, MUX, N, SIGMAX, MUY, M, SIGMAY, 
RLHAT(NINTCREAL( REP)**( 1. 0-ALPHA))), * eHAT(GLOSE), 
REAL( CLOSE) / 1000. 0, ERRBAR, ERRSV 
1900 FORMATC' SIMULATION: yheats "MUX: »F5.1, T16, 'N: | 12, 
T35,'SIGNAX: ' Fy. 1,/,'MUY: "1 F7.3,T16, 'M: so iv ae 
T35,'SIGMAY: ',F4.1, /;' TRUE R: ',F7.5, 
T35,'RLHAT: ' bre 5 i i CLOSEST RLHAT: '\F7.5, 
T35,'TRUE CONFIDENCE LEVEL: ',F5.3,/, 
‘MEAN ERROR: ',F7.5,/, 
"VARIANCE OF ERROR: sF7.5,///) 


++ 


4 
+ 
+ 
4 
+ 
ae 
2000 CONTINUE 


STOP 
END 
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APPENDIX E. 


FILE: SETUNQ DATA 


16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 


16807 93943 © 


X 
MEAN 


ae 
SEED SEED 


iofoto efoLe Lele LokoLoko kokoro 








SIMULATION PARAMETER SETS FOR UNEQUAL 


VARIANCES CASE 
(FOR APPROXIMATE PROCEDURE) 


7.133 1.0 2.01020 .900 1 
7.133 1.0 2.0 25 35 .900 2 
7.133 1.0 2.0 75 50 .900 3 
247.142 10.0 40.0 10 20 .900 4 
247.142 10.0 40.0 25 35 .900 5 
247.142 10.0 40.0 75 50 .900 6 
6.322 1.0 2.010 20 .950 7 
6.322 1.0 2.0 25 35 .950 8 
6.322 1.0 2.075 50 .950 9 
232.175 10.0 40.0 10 20 .950 10 
232.175 10.0 40.0 25 35 .950 ili 
232.175 10.0 40.0 75 50 .950 12 
4,799 1.0 2.01020 .990 13 
4.799 1.0 2.0 25 35 .990 14 
4.799 1.0 2.075 50 .990 15 
204.097 10.0 40.0 10 20 .990 16 
204,097 10.0 40.0 25 35 .990 17 
204.097 10.0 40.0 75 50 .990 18 


Y xX Y 
MEAN STDV STDV N M R 








FILE: SUQ? DATA (FOR COMPARISON SIMULATIONS). 


16807 93943 
16807 93943 
16807 93943 
16807 93943 


16807 93943 $ 
16807 93943 


16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 
16807 93943 


coooocooocooc*oeoooeo°ce 


7. 
7. 
7 


133 
133 
133 


142 
. 142 
. 142 
. 322 
. 322 
. 322 
. 175 
.175 
.175 


10 
70 


Ce ee ee en ee ne ee ee ey 


SEED- SEED 


MEAN 


NEAN 


1.0 2.0 
1.0 2.0 
1.0 2.0 

10.0 40.0 

10.0 40.0 

10.0 40.0 
1.0 2.0 
1.0 2.0 
1.0 2.0 

10.0 40.0 

10.0 40.0 

10.0 40.0 
1053240 
1.0 2.0 
1.0 2.0 

10.0 40.0 

10.0 40.0 

10.0 40.0 

F4.1 F4.1 

x Y 

STDV STDV 





4§ 


- 900 
. 900 
. 900 
. 900 
. 900 
- 900 
. 950 
- 950 
- 950 
. 950 
- 950 
. 950 
.990 
- 990 
. 990 
.990 
. 990 
. 990 
F5.3 
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APPENDIX F. FORTRAN CODE FOR SIMULTANEOUS COMPARISON 
SIMULATION OF APPROXIMATE PROCEDURE VS. 


NONPARAMETRIC PROCEDURE - UNEQUAL VARIANCES 
PROGRAN CONUQ7 


dedededaevevetedede % tedededededodededcavedeakdedededcsededededeakedeacaededeveaedede dedededevedededeavakdeveden de aeveseae saved severe devessge 
% THIS PROGRAM IS TO COMPUTE THE LOWER CONFIDENCE LIMITS OF * 
* R = P (X > Y) SIMULTANEOUSLY FOR BOTH THE APPROXIMATE PRO- ¥ 
* AND THE NONPARAMETRIC PROCEDURE, WHERE X, Y ARE INDEPENDENTLY % 
* NORMALLY DISTRIBUTED WITH UNKNOWN MEANS AND UNKNOWN AND UNEQUAL  * 


%* —-VARTANCES v 
sk ve 
ve ALPHA - NOMINAL CONFIDENCE LEVEL ¥% 
¥e ANORDF - IMSL FUNCTION FOR NORMAL PROBABILITY ve 
* ANORIN - IMSL FUNCTION FOR INVERSE NORMAL CDF %* 
ve ASVX - SAMPLE VARIANCE OF X DEVIDED BY SAMPLE SIZE % 
* ASVY - SAMPLE VARIANCE OF Y DEVIDED BY SAMPLE SIZE oo 
% BIGU - MANN-WHITNEY STATISTIC * 
* CASE - NUMBER OF TESTING ¥ 
* EPS - WIDTH OF THE CONFIDENCE BOUND ¥ 
ve EROR - EROR OF ESTIMATION (APPROXI PROCEDURE) * 

bi ¥ ERORN ~ EROR OF ESTIMATION (NONPARA PROCEDURE) * 
* ERBAR - MEAN OF EROR (APPROXI PROCEDURE) * 
* ERBARN - MEAN OF EROR (NONPARA PROCEDURE) % 
4 * ERSSQ ~ SUN OF SQUARE OF EROR (APPROXI PROCEDURE) * 
¥ ERSSQN - SUN OF SQUARE OF EROR (NONPARA PROCEDURE) * 
ve ERSUM ~ SUM. OF EROR (APPROXI PROCEDURE) * 
* ERSUMN ~ SUM OF EROR (NONPARA PROCEDURE) * 
ERSV ~ SAMPLE VARIANCE OF EROR (APPROXI PROCEDURE) ve 
ve ERSVN - SAMPLE VARIANCE OF EROR (NONPARA PROCEDURE) ve 
¥ CLOSE - INDEX OF THE CLOSEST ESTIMAT (APPROXI PROCEDURE) % 
* CLOSEN - INDEX OF THE CLOSEST ESTIMAT (NONPARA PROCEDURE) % 
ve DELTAH - DELTAHAT; A ESTIMATOR OF DELTA ve 
ve DF ~ DEGREES OF FREEDOM OF DELTAHAT ¥ 
* LNORM - RANDOM NUMBER GENERATOR FOR NORMAL VARIATES ve 
* M - SAMPLE SIZE OF RANDOM VARIABLE X * 
° MUX - POPULATION MEAN OF X % 
* MUY - POPULATION MEAN OF Y ¥ 
* N ~ SAMPLE SIZE OF RANDON VARIABLE Y ve 
* NU - THE SMALLER OF SAMPLE SIZES * 
“ R - REAL RELIABILITY oo 
* RB - CONFIDENCE BOUND OF THE RELIABILITY ve 
* REP ~ REPETITION OF SIMULATIONS * 
& % RLHAT - LOWER CINFIDENCE LIMIT OF RELIABILITY % 
* RTILD ~ POINT ESTIMATOR OF THE RELIABILITY ve 
* SDHAT - SAMPLE STANDARD DEVIATION OF THE DELTAHAT ve 
we SIGMAX - STANDARD DEVIATION OF X * 
, * SIGMAY - STANDARD DEVIATION OF Y % 
* SUNSGX ~ SUM OF SQUARE OF RANDOM SAMPLE OF X ¥e 
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* SUMSQY ~- SUM OF SQUARE OF RANDOM SAMPLE OF Y 
7 SUMX ~ SUN OF RANDOM SAMPLE OF X ui 
* SUNY ~ SUM OF RANDOM SAMPLE OF Y * 
* SVX ~ SAMPLE VARIANCE OF X * 
* SVY - SAMPLE VARIANCE OF Y * 
% T - PERCENTILE OF THE T DISTRIBUTION ae 
* TIN - IMSL FUNCTION TO. COMPUTE PERCENTILE OF T DISTRIBUTION * 
* VARHAT - SAMPLE VARIANCE OF DELTAHAT a 
* XBAR - SAMPLE MEAN OF X * 
: YBAR - SAMPLE MEAN OF Y 7 
x REQUIRED EXTERNAL FUNCTIONS: LNORM, TIN, ANORDF, ANORIN + 
* * 


denvkdckidedeidokkicdadcdcionkioikivdevektelcdkistedcseacdedeledwededescdsesede se de vededke dese dese de devede sere seve ger 357% 
INTEGER REP, -CASE 
REAL ALPHA 
PARAMETER( REP=1000) 
PARAMETER( CASE=18) 
PARAMETER( ALPHA=0. 05) 
INTEGER I, I1, J, U, V, XSEED, YSEED, N, M, CLOSE, 
+ A, B, AI, B1, CLOSEN 
REAL RM, RN, MUX, MUY, XBAR, YBAR, SIGMAX, SIGMAY, R, 
+ X(100), YC 100), X10 100), ¥1(100), SUMSQX, SUNSQY, SUMX, SUMY, 
+ SVX, SVY, ASVX, ASVY, DF1, DF, T, TIN, DELTAH, VARHAT, 
# SDHAT, RLHAT(REP), ANORDF, TEMP, DIFF, EROR, ERSUM, ERSSQ, 
+ ERBAR, ERSV, 
+ BIGU, RTILD, NU, EPS, ANORIN, ERN, ERSUNN, ERSSQN, ERBARN, 
+ ERSVN, DFF, RB(REP) 


CALL EXCMS('FILEDEF 12 DISK SUQ7 DATA Al') 
CALL EXCMS('FILEDEF 18 DISK AUQ7 DATA Al') 
WRITE (18,7) ‘UNEQUAL VARIANCES' 
DO 2000 I=1, CASE 
READ (12,100) XSEED, YSEED, MUX, MUY, SIGMAX, SIGMAY, N, M, R 
100 FORMAT (15,1X,15,1X,F5. 1,1X,F7. 3,1X,F4. 1,1X,F4.1,1X,12,1X,12,1X, 
+F5. 3) 
RN = REAL(N) 
RM = REAL(N) 
ERSUM = 0.0 
ERSSQ = 0.0 
ERSUNN = 0.0 
ERSSQN = 0.0 
DO 1000 J=1, REP 
CALL LNORM(XSEED, X, N, 2, 0) 
CALL LNORM(YSEED, Y, M, 2, 0) 
SUMSQX = 0.0 
SUMSQY = -0.0 
SUMX = 0.0 
SUMY = 0.0 
DO 200 U= 1, N 
X1(U) = X(U) %* SIGMAX + MUX 
SUNSQX= SUMSQX + X1(U) * X1(U) 
SUMX = SUMX + X1(U) 
200 CONTINUE 
XBAR= SUMX / RN 
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DO 300 V=1, M 
Y1(V) = ¥(V) * SIGMAY + MUY 
SUMSQY = SUNSQY + Y1(V) * Y1(V) 
SUMY = SUMY + Y1(V) 


300 CONTINUE 
¥ YBAR = SUMY / RN 
SVX = (SUMSQX - RN * XBAR * XBAR) / (RN - 1.0) 
SVY = (SUMSQY - RM * YBAR * YBAR) / (RM - 1.0) 
+ ASVX = SVX / RN 


ASVY = SVY / RM 
DF1 = (SVX + SVY) * (SVX + SVY) / ( SVX*SVX / (RN-1.0) + 
+ ‘SVY*SVY / (RM-1.0) ) 
DF = ANINT(DF1) 
T = TIN( 1.0-ALPHA, DF) 
DELTAH = (XBAR ~ YBAR) / SQRT(SVX + SVY) 
VARHAT = (ASVX + ASVY) / (SVX + SVY) + 
(XBAR - YBAR) * (XBAR - YBAR) / 
( 2. 0*(SVX + SVY)3 ) % 
( (SVX#"2) / (RN-1.0) + (SVY**2) / (RM-1.0) ) 
SDHAT = SQRT(VARHAT) 
RLHAT(J) = ANORDF(DELTAH - T %* SDHAT) 
EROR = R - RLHAT(J) 
ERSUM = ERSUM + EROR 
ERSSQ = ERSSQ + EROR * EROR 


“ PROCEDURE FOR NONPARAMETRIC 


b+ + 


BIGU = 0.0 
DO 500A =1, N 
‘ DO 400 B = 1, M 
IF (X1(A) .GT. Yi(B)) BIGU = BIGU + 1.0 
400 CONTINUE 
500 CONTINUE 


RTILD = BIGU / (RN * RM) 
NU = MIN (RN, RM) 
EPS = 1.0 / SQRT(4.0 * NU) * ANORIN(1.0 - ALPHA) 
RB(J) = RTILD - EPS 
ERN = R - RB(J) 
ERSUMN = ERSUMN + ERN 
ERSSON = ERSSQN + ERN ** ERN 
1000 CONTINUE 


ERBARN = ERSUMN / REAL(REP) 
ERSVN = ( ERSSQN - REALCREP) * ERBARN * ERBARN ) / REAL(REP-1) 
DFF = 2.0 
DO 1300 Al = 1, REP 
DO 1200 Bl = Al + 1, REP 
IF (RB(B1) . LT. RB(A1)) THEN 
TEMP = RB(A1) 


RB(AL) = RB(B1) 
RB(31) = TEMP 
‘ ENDIF 
1200 CONTINUE 


IF (((R-RB(A1)) .GE. 0.1E-6) . AND. 
+  ((R-RB(A1)) .LE. DFF)) THEN 
: DFF = R - RB(A1) 
CLOSEN = Al 
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ENDIF 
1300 CONTINUE 
v PROCEDURE OF PARAMETRIC 
ERBAR = ERSUM / REALCREP) 
ERSV = ( ERSSQ ~ REAL(REP) * ERBAR * ERBAR ) / REAL(REP-1) 
DIFF = 2.0 
DO 1800 11=1, -REP 
DO 1500 J=1I1+1, REP 
IF (RLHAT(J) . LT. RLHAT(I1)) THEN 
TEMP = RLHAT(1I1) 
RLHAT(I1) = RLHAT(J) 
RLHAT(J) = TEMP 
ENDIF 
1500 CONTINUE 
IF ((C(R-RLHAT(I1)) .GE. 0.1E-6) . AND. 
+  ((R-RLHAT(I1)) .LE. DIFF)) THEN 
DIFF =R - RLHAT( 11) 
CLOSE = I1 
ENDIF 
1800 CONTINUE 
WRITE (18,1500) I, N, M, R, SIGMAX, SIGMAY 
WRITE (18,1910) RLHAT( NINTCREAL(REP)*( 1. 0-ALPHA))), ERBAR 
WRITE (18, 1920) RB(NINT(REAL( REP)*(1. O-ALPHA))), rere 
1900 FORMAT('CASE: ',12,/,'N: ', 12, T16, 'M: ‘, 12, J, 'R: ', P4.3, 
+ Tl6, "STIGMA. x: ', F4.1, T35, ‘SIGMA Y: ',F4.1) 
1910 FORNAT('< PARAMETRIC Bho Ps "RLHAT: ', FS. 4, T16, "ERROR: ', F5.4) 
1920 FORMAT('< NONPARAMETRIC >t? /; 'RLHAT: ' ee 4, T16, "ERROR: an 
+ F5.4,//) 
2000 CONTINUE ; 
STOP 
END 
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